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Introduction 

A matrix is totally positive (resp. totally nonnegative) if all its minors are pos- 
itive (resp. nonnegative) real numbers. The first systematic study of these classes 
of matrices was undertaken in the 1930s by F. R. Gantmacher and M. G. Krein 
|20[ ||^] , who established their remarkable spectral properties (in particular, 
an n X n totally positive matrix x has n distinct positive eigenvalues). Earlier, 
I. J. Schoenberg [SJ discovered the connection between total nonnegativity and 
the following variation-diminishing property: the number of sign changes in a vec- 
tor does not increase upon multiplying by x. 

Total positivity found numerous applications and was studied from many dif- 
ferent angles. An incomplete list includes: oscillations in mechanical systems (the 
original motivation in [ p2[ ) , stochastic processes and approximation theory |2^. , 
Polya frequency sequences |2^, ^] , representation theory of the infinite symmetric 
group and the Edrei-Thoma theorem |l^, planar resistor networks |n|], uni- 
modality and log-concavity Ji^ ] , and theory of immanants ]43| ] . Further references 
can be found in S. Karlin's book [^8| and in the surveys P,]5|^|3^]. 

In this paper, we focus on the following two problems: 

(i) parametrizing all totally nonnegative matrices; 

(ii) testing a matrix for total positivity. 

Our interest in these problems stemmed from a surprising representation-theoretic 
connection between total positivity and canonical bases for quantum groups, dis- 
covered by G. Lusztig Q (cf. also the surveys jn], HU). Among other things, he 
extended the subject by defining totally positive and totally nonnegative elements 
for any reductive group. Further development of these ideas in [|[ [|, [l5[ [lTj aims at 
generalizing the whole body of classical determinantal calculus to any semisimplc 
group. 

As it often happens, putting things in a more general perspective shed new 
light on this classical subject. In the next two sections of this paper, we provide 
self-contained proofs (many of them new) of the fundamental results on problems 
(i)-(ii), due to A. Whitney @, C. Locwncr |§, C. Crycr §, |o), and M. Gasca and 
J. M. Pena |^|. The rest of the paper presents more recent results obtained in Jl5[ : 
a family of efficient total positivity criteria, and explicit formulas for expanding a 
generic matrix into a product of elementary Jacobi matrices. These results and 
their proofs can be generalized to arbitrary semisimple groups [l5|], but we do 
not discuss this here. 

Our approach to the subject relies on two combinatorial constructions. The first 
one is well known: it associates a totally nonnegative matrix to a planar directed 
graph with positively weighted edges (in fact, every totally nonnegative matrix can 
be obtained in this way Q). Our second combinatorial tool was introduced in p!B[ ; 
it is a particular class of colored pseudoline arrangements that we call the double 
wiring diagrams. 
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Planar networks 

To the uninitiated, it might be unclear that totally positive matrices of arbitrary 
order exist at all. As a warm-up, we invite the reader to check that every matrix 
given by 

d dh dhi 

(1) bd bdh + e bdhi + eg + ei 

abd abdh + ae + ce abdhi + (a + c)e(g+i) + f 

where the numbers a, b, c, d, e, /, g, h, i are positive, is totally positive. It will follow 
from the results below that every 3x3 totally positive matrix has this form. 

We will now describe a general procedure that produces totally nonnegative 
matrices. In what follows, a planar network (T, iS) is an acyclic directed planar 
graph r whose edges e are assigned scalar weights uj{e). In all of our examples 
(Figures we assume the edges of T directed left to right. Also, each of our 

networks will have n sources and n sinks, located at the left (resp. right) edge of 
the picture, and numbered bottom-to-top. 




Figure 1. A planar network 



The weight of a directed path in T is defined as the product of the weights of 
its edges. The weight matrix x(T,uj) is an n x n matrix whose (i,j)-entry is the 
sum of weights of all paths from the source i to the sink j. For example, the weight 
matrix of the network in Figure [l] is given by ([!]) . 

The minors of the weight matrix of a planar network have an important com- 
binatorial interpretation, which can be traced to B. Lindstrom J3(^] and further to 
S. Karlin and G. McGregor p9[ | (implicit), and whose many applications were given 
by I. Gessel and G. X. Vicnnot g, 

In what follows, A/ ! j(x) denotes the minor of a matrix x with the row set / and 
the column set J. 

The weight of a collection of directed paths in T is defined to be the product of 
their weights. 

Lemma 1. (Lindstrom's Lemma) A minor A/^j of the weight matrix of a planar 
network is equal to the sum of weights of all collections of vertex-disjoint paths that 
connect the sources labeled by I with the sinks labeled by J . 

To illustrate, consider the matrix x in ([!]). We have, e.g., A23.23(x) = bedegh + 
bdfh + fe, which also equals the sum of the weights of the three vertex-disjoint path 
collections in Figure |l| that connect sources 2 and 3 to sinks 2 and 3. 

Proof. It suffices to prove the lemma for the determinant of the whole weight 
matrix x = x(T,uj), i.e., for the case I = J = [l,nj. Expanding the determinant, 
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we obtain 



(2) 



det(x) = ^2 ^2 s S n ( w ) 



the sum being over all permutations w in the symmetric group S n , and over all col- 
lections of paths tt = (jri, . . . , 7r„) such that m joins the source i with the sink w(i). 
Any collection tt of vertex-disjoint paths is associated with the identity permuta- 
tion; hence oj (tt) appears in (||) with the positive sign. We need to show that all 
other terms in (0) cancel out. Deforming T a bit if necessary, we may assume that 
no two vertices lie on the same vertical line. This makes the following involution on 
the non-vertex-disjoint collections of paths well-defined: take the rightmost point 
of intersection of two paths in tt, and switch the parts of these paths lying to the 
right of this point. This involution preserves the weight of tt, while changing the 
sign of the associated permutation w; the corresponding pairing of terms in (|^) 
provides the desired cancellation. □ 



Corollary 2. If a planar network has nonnegative real weights, then its weight 
matrix is totally nonnegative. 

An aside: note that the weight matrix of the network 

n* * * * • n 




(with unit edge weights) is the "Pascal triangle'' 



which is totally nonnegative by Corollary |2j. Similar arguments can be used to show 
total nonnegativity of various other combinatorial matrices, such as the matrices 
of q-binomial coefficients, Stirling numbers of both kinds, etc. 

We call a planar network T totally connected if, for any two subsets I, J C [1, n] of 
the same cardinality, there exists a collection of vertex-disjoint paths in T connecting 
the sources labeled by I with the sinks labeled by J. 

Corollary 3. If a totally connected planar network has positive weights, then its 
weight matrix is totally positive. 

For any n, let Tq denote the network shown in Figure |^. Direct inspection shows 
that Tq is totally connected. 

Corollary 4. For any choice of positive weights ui(e), the weight matrix x{Tq,uj) 
is totally positive. 
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Figure 2. Planar network T 

It turns out that this construction produces all totally positive matrices; this 
result is essentially equivalent to A. Whitney's Reduction Theorem [Q, and can 
be sharpened as follows. Call an edge of To essential if it either is slanted, or is one 
of the n horizontal edges in the middle of the network. Note that To has exactly n 2 
essential edges. A weighting w of To is essential if u>(e) ^ for any essential edge 
e, and u>(e) = 1 for all other edges. 

Theorem 5. The map u> i— > x(To,u>) restricts to a bijection between the set of all 
essential positive weightings of To and the set of all totally positive n x n matrices. 

The proof of this theorem will use the following notions. A minor A/ t j is called 
solid if both / and J consist of several consecutive indices; if furthermore I U J 
contains 1, then A^j is called initial (see Figure ||). Each matrix entry is the lower- 
right corner of exactly one initial minor; thus the total number of such minors is n 2 . 



Figure 3. Initial minors 

Lemma 6. The n 2 weights of essential edges in an essential weighting u> of To are 
related to the n 2 initial minors of the weight matrix x = x(To,oS) by an invertible 
monomial transformation. Thus an essential weighting u) of Tq is uniquely recovered 
from x. 

Proof. The network Tq has the following easily verified property: for any set / 
of k consecutive indices in [l,n], there is a unique collection of k vertex-disjoint 
paths connecting the sources labeled by [1, k] (resp. by /) with the sinks labeled 
by / (resp. by [1, &]). These paths arc shown by dotted lines in Figure ||, for k = 2 
and / = [3,4]. By Lindstrom's Lemma, every initial minor A of x(Tq,uj) is equal 
to the product of the weights of essential edges covered by this family of paths. 
Notice that among these edges, there is always a unique uppermost essential edge 
e(A) (indicated by the arrow in Figure]^). Furthermore, the map A i— > e(A) is a 
bijection between initial minors and essential edges. It follows that the weight of 
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each essential edge e = e(A) is equal to A times a Laurent monomial in some initial 
minors A' whose associated edges e(A') are located below e. □ 

To illustrate Lemma ^, consider the special case n = 3. The network Tq is 
shown in Figure |l]; its essential edges have the weights a,b,...,i. The weight 
matrix x(To,u>) is given in ([!]). Its initial minors are given by the monomials 

Ai .1 = d Ai : 2 = dh A± t 3 = dhi 
A 2 ,i = bd A12.12 = de Ai 2 ,23 = degh 
A 3> i = abd A23.12 = bade Ai 23 ,i23 = def_ 

where for each minor A, the "leading entry" w(e(A)) is underlined. 

To complete the proof of Theorem ||, it remains to show that every totally 
positive matrix x has the form x(Tq, uj) for some essential positive weighting u). By 
Lemma ^, such an ui can be chosen so that x and x{Tq, lu) will have the same initial 
minors. Thus our claim will follow from the lemma below. 

Lemma 7. A square matrix x is uniquely determined by its initial minors provided 
all these minors are nonzero. 

Proof. Let us show that each matrix entry Xij of x is uniquely determined by the 
initial minors. If i — l or j = l, there is nothing to prove, since is itself an initial 
minor. Assume that min(z, j) > 1. Let A be the initial minor whose last row is i 
and last column is j, and let A' be the initial minor obtained from A by deleting 
this row and this column. Then A = A'xij + P, where P is a polynomial in the 
matrix entries x^ji with ^ and i' < i, j' < j. Using induction on i + j, 

we can assume that each xyji that occurs in P is uniquely determined by the initial 
minors, so the same is true for x^j — (A — P)/A'. This completes the proofs of 
Lemma Q and Theorem ||. □ 

Theorem ^| describes a parametrization of totally positive matrices by n 2 -tuples 
of positive reals, providing a partial answer (one of the many possible, as we will 
see) to the first problem stated in the introduction. The second problem — that of 
testing total positivity of a matrix — can also be solved using this theorem, as we 
will now explain. 

An n x n matrix has altogether ( 2 ? ™) — 1 minors. This makes it impractical to 
test positivity of every single minor. It is desirable to find efficient criteria for total 
positivity that would only check a small fraction of all minors. 

a b 
c d 

entries and the determinant A = ad — be. To test that x is totally positive, it is 
enough to check positivity of a, b, c, and A; then d = (A + be) /a > 0. 

The following theorem generalizes this example to matrices of arbitrary size; it 
is a direct corollary of Theorem || and Lemmas |(3| and 0. 

Theorem 9. A square matrix is totally positive if and only if all its initial minors 
(see Figure are positive. 

This criterion involves n 2 minors, and it can be shown that this number cannot 
be lessened. Theorem || was proved by M. Gasca and J. M. Peha p3|, T heorem 4.1] 
(for rectangular matrices); it also follows from C. Cryer's results in |9|. Theorem 
is an enhancement of the 1912 criterion by M. Fekete H], who proved that the 
positivity of all solid minors of a matrix implies its total positivity. 



Example 8. A 2 x 2 matrix x = 



has ( 2 ) — 1 = 5 minors: four matrix 
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Theorems of Whitney and Loewner 



In this paper, we shall only consider invertible totally nonnegative nxn matrices. 
Although these matrices have real entries, it is convenient to view them as elements 
of the general linear group G = GL n (C). We denote by G>o (resp. G>o) the set of 
all totally nonnegative (resp. totally positive) matrices in G. The structural theory 
of these matrices begins with the following basic observation, which is an immediate 
corollary of the Binet-Cauchy formula. 

Proposition 10. Both G>o and G>o are closed under matrix multiplication. Fur- 
thermore, if x £ G>o and y S G>o then both xy and yx belong to G>o ■ 

Combining this proposition with the foregoing results, we will prove the following 
theorem of A. Whitney p6[ |. 

Theorem 11. (Whitney's Theorem) Every invertible totally nonnegative matrix is 
the limit of a sequence of totally positive matrices. 

Thus G>o is the closure of G>o in G. (The condition of invertibility in Theo- 
rem [n] can in fact be lifted.) 

Proof. First let us show that the identity matrix / lies in the closure of G>o ■ 
By Corollary ^, it suffices to show that / = limjv^oo x(Tq,ujn) for some sequence 
of positive weightings lon of the network Tq . Note that the map ui i— ► x(Tq,u>) is 
continuous, and choose any sequence of positive weightings that converges to the 
weighting luq defined by uJo(e) = 1 (resp. 0) for all horizontal (resp. slanted) edges e. 
Clearly, x^q^loq) = /, as desired. 

To complete the proof, write any matrix x e G>o as x = lini/v— >oo x ■ x(Tq,lon) , 
and notice that all matrices x ■ x(To,wn) are totally positive by Proposition |l0[ □ 

The following description of the multiplicative monoid G>o was first given by 
C. Loewner 32 under the name "Whitney's Theorem"; it can indeed be deduced 
from pi. 



Theorem 12. (Loewner- Whitney Theorem) Any invertible totally nonnegative 
matrix is a product of elementary Jacobi matrices with nonnegative matrix entries. 

Here "elementary Jacobi matrix" is a matrix xGG that differs from / in a single 
entry located either on the main diagonal, or immediately above or below it. 
Proof. We start with an inventory of elementary Jacobi matrices. Let Eij denote 
the nxn matrix whose (i, j)-entry is 1 while all other entries are 0. For t 6 C and 
i = 1 , . . . , n — 1 , let 



Xi(t) = I + tE lA+1 











and 



x- i (t)=I + tE i+1 , i = (x t (t)y 
(the transpose of Xi{t)). Also, for i = l,...,n and t ^ 0, let 

x @ (t) = I + (t - l)E iti , 
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the diagonal matrix with ith diagonal entry equal to t and all other diagonal entries 
equal to 1. Thus elementary Jacobi matrices are precisely the matrices of the form 
Xi(t), xj(t), and x©(t). An easy check shows that they are totally nonnegative for 
any t > 0. 

For any word i = (ii, . . . , i{) in the alphabet 

(3) A={l,...,n-1, (£,..., ® ,T,...,ir=T} 
we define the product map Xi : (C \ {0})' — > G by 

(4) Xi(ti, ...,ti) = x tl (ti) ■ ■ -Xi^U) . 

(Actually, Xi(t\, . . . ,ti) is well defined as long as the right-hand side of (^) does not 
involve any factors of the form x©(0).) To illustrate, the word i = (T) 1(2)1 gives 



rise to 



Xi(ti,t2,t 3 ,t 4 ) = 



h 
o 

h 

t-2 





1 




1 









*2 1 







h 





u 
1 



We will interpret each matrix a!i(ti, . . . , t{) as the weight matrix of a planar network. 
First note that any elementary Jacobi matrix is the weight matrix of a "chip" of 
one of the three kinds shown in Figure [I| In each "chip," all edges but one have 
weight 1; the distinguished edge has weight t. Slanted edges connect horizontal 
levels i and i + 1, counting from the bottom; in all examples in Figure 0, % = 2. 





Xi(t) Xj(t) x@(t) 

Figure 4. Elementary "chips" 

The weighted planar network (r(i), w{t\, . . . , ij)) is then constructed by con- 
catenating the "chips" corresponding to consecutive factors x% h {tk), as shown in 
Figure [|. It is easy to see that concatenation of planar networks corresponds to 
multiplying their weight matrices. We conclude that the product Xi(t\, . . . ,ti) of 
elementary Jacobi matrices equals the weight matrix x(T(i), oj{t\, . . . , £;)). 




i = 21@21®21(2) 
Figure 5. Planar network T(i) 
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In particular, the network (Tq,lj) appearing in Figure |^ and Theorem |5| (more 
precisely, its equivalent deformation) corresponds to some special word i max of 
length n 2 ; instead of defining i max formally, we just write it for n = 4: 

W = (3, 2, 3, T, 2, 3,®, ©, ©, @, 3, 2, 3, 1, 2, 3) . 

In view of this, Theorem || can be reformulated as follows. 

Theorem 13. The product map Xi max restricts to a bijection between n 2 -tuples of 
positive real numbers and totally positive n x n matrices. 



We will prove the following refinement of Theorem 12, which is a reformulation 
of its original version j}2| . 

Theorem 14. Every matrix x E G>o can be written as x = £i max (£i, ■ ■ ■ , £ r i 2 )> f or 
some ti, . . . , t„2 > 0. 

(Since x is invertible, we must in fact have tk > for n ("~ 1 ) < k < 21iZLhil ] i. e-) 
for those indices fc for which the corresponding entry of i ma x is of the form ©.) 
Proof. The following key lemma is due to C. Cryer ||. 

Lemma 15. The leading principal minors An y u w of a matrix x E G>o are 
positive, for fe = l, . . . , n. 

Proof. Using induction on n, it suffices to show that A\i t n—i],[i,n-i\{ x ) > 0- Let 
A % '^{x) (resp. A M JJ (x)) denote the minor of x obtained by deleting the row i 
and the column j (resp. rows i and i', and columns j and j'). Then, for any 
1 < i < i' < n and 1 < j < j' < n, one has 

(5) A i '' j '(x)A i ' j (a?) - A J ' J (x)A ij '(x) = det(x) A u '^\x) 

as an immediate consequence of Jacobi's formula for minors of the inverse ma- 
trix (see, e.g., 0, Lemma 9.2.10]). The determinantal identity (||) was proved by 
P. Desnanot as early as in 1819 (see (3?j, pp. 140-142]); it is sometimes called "Lewis 
Carroll's identity," due to the role it plays in C. L. Dodgson's condensation method 
(IH pp. 170-180]. 

Now suppose that A n ' n (x) = for some x E G>o ■ Since x is invertible, we have 
A hn (x) > and A n ' j (x) > for some indices i,j < n. Using (||) with i' = j' = n, 
we arrive at a desired contradiction by 

> -A n ' j (x)A i ' n (x) = det(x)A in ' jn (x) > . □ 

We are now ready to complete the proof of Theorem |lj. Any matrix x E G>o 
is by Theorem ^ a limit of totally positive matrices x jv , each of which can by 
Theorem |l3| be factored as xn = x i m ^(t[ N ^ ■ ■ ■ > with all t^' positive. It 
suffices to show that the sequence sjv = Sfc=i*fe converges; then the standard 
compactness argument will imply that the sequence of vectors (tj , . . . , t^P) con- 
tains a converging subsequence, whose limit (ti,. . ■ ,t n z) will provide the desired 
factorization x = 2ii max (ti, ■ ■ ■ , ^n 2 )- To see that (sjv) converges, we use the explicit 
formula 

^[l,i],[l,i](^jv) 



S N 



> A[i ) i_i]u{i+i},[i,»](icjv) + A[ lji ] j [ lj j_ 1 ] U { i+1 }(a;jv) 



+ A [M]j[M] (xjv) 
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(to prove this, compute the minors on the right with the help of Lindstrom's Lemma 
and simplify). Thus sn is expressed as a Laurent polynomial in the minors of xn 
whose denominators only involve leading principal minors A[i w n,fc] . By Lemma [l5| 
as x n converges to x, this Laurent polynomial converges to its value at x. This 
completes the proofs of Theorems |lj and |l4|. □ 

Double wiring diagrams and 
total positivity criteria 

We will now give another proof of Theorem [)[ which will include it into a family of 
"optimal" total positivity criteria that correspond to combinatorial objects called 
double wiring diagrams. This notion is best explained by an example, such as 
the one given in Figure g. A double wiring diagram consists of two families of n 
piecewise-straight lines (each family colored with one of the two colors), the crucial 
requirement being that each pair of lines of like color intersect exactly once. 




Figure 6. Double wiring diagram 

The lines in a double wiring diagram are numbered separately within each color. 
We then assign to every chamber of a diagram a pair of subsets of the set [1, n] = 
{1, . . . , n}: each subset indicates which lines of the corresponding color pass below 
that chamber; see Figure 



123,123 




Figure 7. Chamber minors 

Thus every chamber is naturally associated with a minor A/.,/ of an n x n 
matrix x = (xij) (we call it a chamber minor) that occupies the rows and columns 
specified by the sets I and J written into that chamber. In our running example, 
there are 9 chamber minors (the total number is always n 2 ), namely X31 , X32 , £12 , 
X13 , A 2 3,i2 , Ai 3 ,i2 , Ai 3j 23 , Ai 2 , 2 3 , and Ai23,i2 3 = det(x). 

Theorem 16. Jl5[ Every double wiring diagram gives rise to the following crite- 
rion: an n xn matrix is totally positive if and only if all its n 2 chamber minors are 
positive. 

The criterion in Theorem [9] is a special case of Theorem [lfj], and arises from the 
"lexicographically minimal" double wiring diagram, shown in Figure || for n = 3. 
Proof. We will actually prove the following statement that implies Theorem 061 
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123,123 




Figure 8. Lexicographically minimal diagram 



Theorem 17. Jl5[ Every minor of a generic square matrix can be written as a 
rational expression in the chamber minors of a given double wiring diagram, and 
moreover this rational expression is subtraction- free, i.e., all coefficients in the 
numerator and denominator are positive. 

Two double wiring diagrams are called isotopic if they have the same collections 
of chamber minors. The terminology suggests what is really going on here: two 
isotopic diagrams have the same "topology." From now on, we will treat such 
diagrams as indistinguishable from each other. 

We will deduce Theorem ^ from the following fact: any two double wiring 
diagrams can be transformed into each other by a sequence of local "moves" of 
three different kinds, shown in Figure ^. (This is a direct corollary of a theorem of 
G. Ringcl |3{| . It can also be derived from the Tits theorem on reduced words in 
the symmetric group; cf. (^) — (j8|) below.) 





B B 

D D 

Figure 9. Local "moves" 

Notice that each local move exchanges a single chamber minor Y with another 
chamber minor Z, and keeps all other chamber minors in place. 

Lemma 18. Whenever two double wiring diagrams differ by a single local move 
of one of the three types shown in Figure Q the chamber minors appearing there 
satisfy the identity AC + BD = YZ . 

The three-term determinantal identities of Lemma [l8] are well known, although 
not in this disguised form. The last of these identities is nothing but the identity 
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, applied to various submatrices of an n x n matrix. The identities corresponding 
to the top two "moves" in Figure ^| are special instances of the classical Grassmann- 
Pliicker relations (see, e.g., (lj, (15.53)]), and were obtained by P. Desnanot along- 
side Q in the same 1819 publication that we mentioned before. 

Theorem [l^ is now proved as follows. We first note that any minor appears 
as chamber minor in some double wiring diagram. Therefore, it suffices to show 
that the chamber minors of one diagram can be writen as subtraction-free rational 
expressions in the chamber minors of any other diagram. This is a direct corollary of 
Lemma |l8] combined with the fact that any two diagrams are related by a sequence 
of local moves: indeed, each local move replaces Y by (AC + BD)/Z, or Z by 
(AC + BD)/Y. □ 

Implicit in the above proof is an important combinatorial structure lying behind 
Theorems [l6] and ^ the graph <f>„ whose vertices are the (isotopy classes of) double 
wiring diagrams, and whose edges correspond to local moves. The study of is an 
interesting problem in itself. The first nontrivial example is the graph $3 shown in 
Figure [h]. It has 34 vertices, corresponding to 34 different total positivity criteria. 
Each of these criteria tests 9 minors of a 3x3 matrix. Five of these minors, viz. X31 , 
CC13 , A23,i2 , Ai2 ; 23 j and det(x), correspond to the "unbounded" chambers that lie 
on the periphery of every double wiring diagram; they are common to all 34 criteria. 
The other four minors correspond to the bounded chambers, and depend on the 
choice of a diagram. For example, the criterion derived from Figure [jj involves 
"bounded" chamber minors A3 2 , Ai^ , A1312 , and A13.23 ■ In Figure [II], each 
vertex of $3 is labeled by the quadruple of "bounded" minors that appear in the 
corresponding total positivity criterion. 



We suggest the following refinement of Theorem 17 



Conjecture 19. Every minor of a generic square matrix can be written as a Lau- 
rent polynomial with nonnegative integer coefficients in the chamber minors of an 
arbitrary double wiring diagram. 

Perhaps more important than proving this conjecture would be to give explicit 
combinatorial expressions for the Laurent polynomials in question. We note a 
case where the conjecture is true, and the desired expressions can be given: the 
"lexicographically minimal" double wiring diagram whose chamber minors are the 
initial minors. Indeed, a generic matrix x can be uniquely written as the product 
Xi max (ti , . . . , t„2 ) of elementary Jacobi matrices (cf. Theorem |l3|) ; then each minor 
of x can be written as a polynomial in the tk with nonnegative integer coefficients 
(with the help of Lindstrom's Lemma), while each tk is a Laurent monomial in the 
initial minors of x, by Lemma [| 

It is proved in ||l5| , Theorem 1.13] that every minor can be written as a Laurent 
polynomial with integer (possibly negative) coefficients in the chamber minors of a 
given diagram. Note, however, that this result, combined with Theorem]!?], does not 
imply Conjecture [l9|, because there do exist subtraction-free rational expressions 
that are Laurent polynomials although not with nonnegative coefficients (e.g., think 
of (p 3 + q 3 )/(p + q)=p 2 - pq + q 2 ). 

The following special case of Conjecture |l9| can be derived from Thm. 3.7.4]. 

Theorem 20. Conjecture holds for all wiring diagrams in which all intersec- 
tions of one color precede the intersections of another color. 
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We do not know an elementary proof of this result; the proof in |3| depends on 
the theory of canonical bases for quantum general linear groups. 
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Digression: Somos sequences. The three-term relation AC + BD = YZ is sur- 
rounded by some magic that eludes our comprehension. We cannot resist mention- 
ing the related problem involving the Somos-5 sequences ]l9| . (We thank Richard 
Stanley for telling us about them.) These arc the sequences 01, a-x, ■ ■ ■ in which any 
6 consecutive terms satisfy this relation: 

(6) 

Each term of a Somos-5 sequence is obviously a subtraction-free rational expression 
in the first 5 terms ai, . . . , as . More can be shown by extending the arguments 
in (1^, ^35|: a„ is actually a Laurent polynomial in a\, . . . , as . This property is truly 
remarkable, given the nature of the recurrence, and the fact that, as n grows, these 
Laurent polynomials become huge sums of monomials involving large coefficients; 
still, each of these sums cancels out from the denominator of the recurrence relation 

»n+5 = (an+lflri+4 + a-n+'l^n+z) / 0>n ■ 

We suggest the following analogue of Conjecture [l9| 

Conjecture 21. Every term of a Somos-5 sequence is a Laurent polynomial with 
nonnegative integer coefficients in the first 5 terms of the sequence. 



Factorization schemes 

According to Theorem [l6], every double wiring diagram gives rise to an "optimal" 
total positivity criterion. We will now show that double wiring diagrams can be 
used to obtain a family of bijective paramctrizations of the set G>o of all totally 
positive matrices; this family will include as a special case the paramctrization in 
Theorem [D| 

We encode a double wiring diagram by the word of length n(n—l) in the alphabet 
{1, . . . ,n — 1, 1, ... ,n — 1} obtained by recording the heights of intersections of 
pseudolines of like color (traced left to right). For example, the diagram in Figure || 
is encoded by the word 2 12 12 1. 

The words that encode double wiring diagrams have an alternative description 
in terms of reduced expressions in the symmetric group S n . Recall that by a 
famous theorem of E. H. Moore [^6| , S n is a Coxeter group of type A n _±, i.e., it is 
generated by the involutions s±, . . . , s„_i (adjacent transpositions) subject to the 
relations SiSj = SjSi for \i — j\ > 2, and SiSjSi = SjSiSi for \i — j\ = 1. A reduced 
word for a permutation w G S n is a word j = (ji, ■ ■ ■ , ji) of the shortest possible 
length I = £(w) that satisfies w = Sj t ■ • ■ Sj l . The number £(w) is called the length 
of w (it is the number of inversions in w). The group S n has a unique element w 
of maximal length: the order-reversing permutation of 1, ... , Ti- 
lt is straightforward to verify that the encodings of double wiring diagrams are 
precisely the shuffles of two reduced words for w , in the barred and unbarred entries 
respectively; cquivalcntly, these are the reduced words for the element (w ,w ) of 
the Coxeter group S n x S n . 

Definition 22. A word i in the alphabet A (see (||)) is called a factorization scheme 
if it contains each circled entry exactly once, and the remaining entries encode 
the heights of intersections in a double wiring diagram. 

Equivalcntly, a factorization scheme i is a shuffle of two reduced words for w (one 
barred and one unbarred) and an arbitrary permutation of the entries ® , . . . , © . 
In particular, i consists of n 2 entries. 
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To illustrate, the word i = 21(3)21(T)21@ appearing in Figure |^ is a factor- 
ization scheme. 

An important example of a factorization scheme is the word i ma x introduced in 
Theorem |l3[ Thus the following result generalizes Theorem [l3|. 

Theorem 23. JT5| For an arbitrary factorization scheme i = (ii, . . . , i n ^), the 
product map x\ given by ^ restricts to a bijection between n 2 -tuples of positive 
real numbers and totally positive n x n matrices. 

Proof. We have already stated that any two double wiring diagrams are connected 
by a succession of the local "moves" shown in Figure ^|. In the language of factor- 
ization schemes, this translates into any two factorization schemes being connected 
by a sequence of local transformations of the form 



(7) 

or of the form 
(8) 



■ ■ ■ l J l ■ ■ ■ 

■ ■ - i j i ■ ■ ■ 



•JiJ--- 
■jij--- 



|i-il = l. 

K-j'I = i, 



• • • ab ■ ■ ■ ■ ■ - ba- ■ 



where (a, b) is any pair of symbols in A different from (i, i ± 1) or (i, i ± 1). (This 
statement is a special case of Tits' theorem [fl5| , for the Coxeter group S n x S n x 
(«5 2 )".) 

In view of Theorem [l3[ it suffices to show that if Theorem ^] holds for some 
factorization scheme i, then it also holds for any factorization scheme i' obtained 
from i by one of the transformations (0)-(P|)- To see this, it is enough to demonstrate 
that the collections of parameters {tk} and {t' k } in the equality 

x h (ii) • • • Xi n2 (t n 2) = x Vi (t[) ■ ■ ■ x^ 2 {t' n2 ) 

are related to each other by (invertible) subtraction-free rational transformations. 
The latter is a direct consequence of the commutation relations between elementary 
Jacobi matrices, which can be found in |l5|, Section 2.2 and (4.17)]. The most 
important of these relations are the following. 
First, for i = 1, . . . , n — 1 and j = i + 1, we have 



X l {tl)x & {t 2 ) XQ)(t 3 ) Xj(ti) = x-(t[) x & {t' 2 ) x@(4) Xi(t' A ) , 



where 



t\t^ti 



The proof of this relation (which is the only nontrivial relation associated with (g|)) 
amounts to verifying that 



1 h 
1 





h o 




i 







1 




t' 2 






o h 




u 


1 




t[ 1 




t' 3 





i A 
o 1 



Also, for any i and j such that = 1, we have the following relation associated 
with (0): 

Xi(tl)Xj(t2)Xi(t 3 ) = Xjit'^Xiit^Xjit^) , 



xjit^xj^xjih) = x-it'^xjit^xjit'a) , 
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T = h + t 3 



where 

— , t' 2 =T, 4 = ^ 

One sees that in the commutation relations above, the formulas expressing the t' k 
in terms of the U are indeed subtraction-free. □ 



Theorem ^3] suggests an alternative approach to total positivity criteria via the 
following factorization problem: for a given factorization scheme i, find the gener- 
icity conditions on a matrix x assuring that x can be factored as 

(9) x = xi(ti, . . . ,t n s) = Xi^ti) ■ ■ -Xi n2 (t n 2) , 

and compute explicitly the factorization parameters tk as functions of x. Then 
the total positivity of x will be equivalent to the positivity of all these functions. 
Note that the criterion in Theorem ^ was essentially obtained in this way: for the 
factorization scheme i max , the factorization parameters t k are Laurent monomials 
in the initial minors of x (cf . Lemma |^) . 

A complete solution of the factorization problem for an arbitrary factorization 
scheme was given in |l5|, Theorems 1.9 and 4.9]. An interesting (and unexpected) 
feature of this solution is that in general, the tk are not Laurent monomials in 
the minors of x; the word i max is quite exceptional in this respect. It turns out, 
however, that the tk are Laurent monomials in the minors of another matrix x' 
obtained from x by the following birational transformation: 

(10) x' = [x T w ] + w (x T y 1 w [w x T ]- . 

Here x T denotes the transpose of x, and w a is the permutation matrix with l's on the 
antidiagonal; finally, y = [y]—[y]o[y]+ denotes the Gaussian (LDU) decomposition 
of a square matrix y provided such a decomposition exists. 
In the special cases n = 2 and n = 3, the transformation x 



x' is given by 



Xi\X 12 x 2 i 



[ 21 



x 2 2 det(x) 



and 



Xll 


Ai2,l3 


1 


X31 Xi 3 


#31 Ai2,23 


X31 




£33Ai2,i2 - det{x) 


X32 


Xl3 A2342 


A23,12 Ai2,23 


A23,12 


1 


£23 


A23,23 


X13 


Ai2,23 


det(x) 



The following theorem provides an alternative explanation for the family of total 
positivity criteria in Theorem [l6|. 

Theorem 24. Jl5[ The right-hand side of (10) is well defined for any x G G>o ; 

moreover, the "twist map" x 1— > x' restricts to a bijection of G>o with itself. 

Let x be a totally positive n x n matrix, and i a factorization scheme. Then 
the parameters ti , . . . , t n 2 appearing in are related by an invertible monomial 
transformation to the n 2 chamber minors (for the double wiring diagram associated 
with 1) of the twisted matrix x' given by ( p^| ) . 

In |l5| , we explicitly describe the monomial transformation in Theorem ^4|, as 
well as its inverse, in terms of the combinatorics of the double wiring diagram. 
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Double Bruhat cells 

Our presentation in this section will be a bit sketchy; details can be found in Jl5[ . 

Theorem ^3] provides a family of bijective (and biregular) parametrizations of 
the totally positive variety G>o by n 2 -tuples of positive real numbers. The totally 
nonncgative variety G>o is much more complicated (note that the map in Theo- 
rem 14 is surjective but not injective). In this section, we show that G>o splits 



naturally into "simple pieces" corresponding to pairs of permutations from S n . 

Theorem 25. Jl5[ Let x 6 G>o be a totally nonnegative matrix. Suppose that a 
word i in the alphabet A is such that x can be factored as x = x\ (t\ , . . . , t m ) with pos- 
itive t\ , . . . , t m , and i has the smallest number of uncircled entries among all words 
with this property. Then the subword of i formed by entries from {1, . . . , n — 1} 
(resp. from {l,...,n— 1}) is a reduced word for some permutation u (resp. v) 
in S n . Furthermore, the pair (u,v) is uniquely determined by x, i.e., does not 
depend on the choice of i. 



In the situation of Theorem |2J], we say that x is of type (u, v). Let G^ C G>o 
denote the subset of all totally nonnegative matrices of type (u, v); thus G>o is the 
disjoint union of these subsets. 

Every subvaricty G"'g has a family of parametrizations similar to those in The- 
orem Generalizing Definition E2], let us call a word i in the alphabet A a 
factorization scheme of type (u, v) if it contains each circled entry exactly once, 
and the barred (resp. unbarred) entries of i form a reduced word for u (resp. v); in 
particular, i is of length £(u) + £(v) + n. 

Theorem 26. Jl5[ For an arbitrary factorization scheme i of type (u,v), the prod- 
uct map x\ restricts to a bisection between (£(u) + £(v) + n)-tuples of positive real 
numbers and totally nonnegative matrices of type (u,v). 

Comparing Theorems and ^3|, we see that 

(11) G>Q w ° = G> , 

i.e., the totally positive matrices are exactly the totally nonnegative matrices of 
type (w ,w ). 

We now show that the splitting of G>o into the union of varieties G^'q is 
closely related to the well-known Bruhat decompositions of the general linear group 
G = GL n . Let B (resp. £?_) denote the subgroup of upper-triangular (resp. lower- 
triangular) matrices in G. Recall (see, e.g., jl|, §4]) that each of the double coset 
spaces B\G/B and BJ\G/ B- has cardinality nl, and one can choose the permuta- 
tion matrices w 6 S n as their common representatives. To every two permutations 
u and v we associate the double Bruhat cell G u,v = BuB n B_vB_\ thus G is the 
disjoint union of the double Bruhat cells. 

Each set G u > v can be described by equations and inequalities of the form A (a;) = 

and/or A(x) ^ 0, for some collection of minors A. (See |l5|, Proposition 4.1] 
or [Q.) In particular, the open double Bruhat cell G w °' w ° is given by non- 
vanishing of all "antiprincipal" minors A[i,i]j n -i+i,n\( x ) an d A[n-t+i,n],[i,i] (x) for 

1 = 1 , . . . , n — 1 . 

Theorem 27. Jl5| A totally nonnegative matrix is of type (u, v) if and only if it 
belongs to the double Bruhat cell G u,v . 
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In view of ([0]), Theorem [27| provides the following simple test for total positivity 
of a totally nonnegative matrix. 

Corollary 28. |23fi A totally nonnegative matrix x is totally positive if and only if 
&li,i],[n-i+i,n](x) and A [ „_ i+lj „ ]j [ lji] (x) ^ for i = 1, . . . , n. 

The results obtained above for G™° Q ' W ° = G>o (as well as their proofs) extend to 
the variety G"'q for an arbitrary pair of permutations u, v £ S n . In particular, the 
factorization schemes for (w, v) (or rather their uncircled parts) can be visualised 
by double wiring diagrams of type (it, v) in the same way as before, except now any 
two pscudolincs intersect at most once, and the lines are permuted "according to u 
and v." Every such diagram has £(u)+£(v) + n chamber minors, and their positivity 
provides a criterion for a matrix x £ G u ' v to belong to G" q" . The factorization 
problem and its solution provided by Theorem extend to any double Bruhat 
cell, with an appropriate modification of the twist map x <— > x' . The details can be 
found in ]l5[ . 

If the double Bruhat cell containing a matrix x £ G is not specified, then testing 
x for total nonnegativity becomes a much harder problem; in fact, every known 
criterion involves exponentially many (in n) minors. (See Q for related complexity 
results.) The following corollary of a result by Cryer Q was given by Gasca and 
Pena @. 

Theorem 29. An invertible square matrix is totally nonnegative if and only if all 
its minors occupying several initial rows or several initial columns are nonnegative, 
and all its leading principal minors are positive. 

This criterion involves 2 n+1 — n — 2 minors, which is roughly the square root of 
the total number of minors. We do not know whether this criterion is optimal. 

Oscillatory matrices 

We conclude the paper by discussing the intermediate class of oscillatory matri- 
ces that was introduced and intensively studied by Gantmacher and Krein ]2(], . 
A matrix is oscillatory if it is totally nonnegative while some power of it is to- 
tally positive; thus the set of oscillatory matrices contains G>o and is contained 
in G>o • The following theorem provides several equivalent characterizations of os- 
cillatory matrices; the equivalence of (a)-(c) was proved in |p2| , while the rest of 
the conditions were given in ||l7f . 

Theorem 30. [^2| , flT] For an invertible totally nonnegative n x n matrix x, the 
following are equivalent: 

(a) x is oscillatory; 

(b) Xi.i+i > and Xi+n > for i = 1, . . . , n— 1; 

(c) x n ~ l is totally positive; 

(d) x is not block-triangular ( cf. Figure [Zj 
(c) x can be factored as x — Xi(b\, . . . , t{), for positive t\, . . . ,ti and a word i that 

contains every symbol of the form i or i at least once; 
(f) x lies in a double Bruhat cell G u ' v , where both u and v do not fix any set 
{1, ...,i},fori = l,.. 
Proof. Obviously, (c) (a) (d). Let us prove the equivalence of (b), (d), 
and (e). By Theorem [l2|, x can be represented as the weight matrix of some planar 
network T(i) with positive edge weights. Then (b) means that sink i + 1 (resp. i) 
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Figure 1 1 . Block-triangular matrices 



can be reached from source i (resp. i + 1), for all i; (d) means that for any i, at 
least one sink j > i is reachable from a source h < i, and at least one sink h < i 
is reachable from a source j > i; and (e) means that T(i) contains positively- and 
negatively-sloped edges connecting any two consecutive levels i and i + 1. These 
three statements are easily seen to be equivalent. 

By Theorem pTj, (c ) <^=^ (f). It remains to show that (c) (c). In view of 
Theorem ^ and (|ll|), this can be restated as follows: given any permutation j 
of the entries 1, . . . , n— 1, prove that the concatenation j™ -1 of n— 1 copies of j 
contains a reduced word for w . Let j' denote the subsequence of j™ -1 constructed 
as follows. First, j' contains all n— 1 entries of j™" 1 which are equal to n—1. Second, 
j' contains the n — 1 entries equal to n-2 which interlace the n — 1 entries chosen 
at the previous step. We then include n — 3 interlacing entries equal to n — 3, etc. 
The resulting word j' of length (™) will be a reduced word for w , for it will be 
equivalent, under the transformations (|J), to the lexicographically maximal reduced 
word j max = (n— 1, n— 2, n— 1, n— 3, n— 2, n— 1, . . .). □ 
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